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ABSTRACT 

The type-I X-ray bursting low mass X-ray binary KS 173 1-260 was recently detected for the first time in quies- 
cence by Wijnands et al., following a r 0Ut b U rst~13 yr outburst which ended in Feb 2001. We show that the emission 
area radius for a H atmosphere spectrum (possibly with a hard power-law component that dominates the emis- 
sion above 3.5 keV) is consistent with that observed from other quiescent neutron star transients, R O0 =23t.f^(d/S 
kpc) km, and examine possible IR counterparts for KS 1731-260. Unlike all other known transient neutron stars 
(NS), the duration of this recent (and the only observed) outburst is as long as the thermal diffusion time of the 
crust. The large amount of heat deposited by reactions in the crust will have heated the crust to temperatures 
much higher than the equilibrium core temperature. As a result, the thermal luminosity currently observed from 
the neutron star is dominated not by the core, but by the crust. This scenario implies that the mean outburst re- 
currence timescale found by Wijnands et al. (~ 200 yr) is a lower limit. Moreover, because the thermal emission 
is dominated by the hot crust, the level and the time evolution of quiescent luminosity is determined mostly by 
the amount of heat deposited in the crust during the most recent outburst (for which reasonable constraints on the 
mass accretion rate exist), and is only weakly sensitive to the core temperature. Using estimates of the outburst 
mass accretion rate, our calculations of the quiescent flux immediately following the end of the outburst agree 
with the observed quiescent flux to within a factor of a few. In this paper, we present simulations of the evolution 
of the quiescent lightcurve for different scenarios of the crust microphysics, and demonstrate that monitoring ob- 
servations (with currently flying instruments) spanning from 1-30 yr can measure the crust cooling timescale and 
the total amount of heat stored in the crust. These quantities have not been directly measured for any neutron star. 
This makes KS 1731-260 a unique laboratory for studying the thermal properties of the crust by monitoring the 
luminosity over the next few years to decades. 

Subject headings: stars: atmospheres — stars: individual (KS 1731-260) — stars: neutron — x-rays: binaries 



1. INTRODUCTION 



Brown, Bildsten & Rutledge (1998, BBR98 hereafter) ar- 



gued that the core of a transiently accreting neutron star (NS) 
(for reviews of transi ent neutron stars, see |Chen et al. 1997 



J> , Campana et al. 1998) is heated by nuclear reactions deep in 



the crust during the accretion outbursts. The cor e is heated to 
a steady state in ~10 4 yr (see |Colpi et al. 2000 for a detailed 
calculation), after which the NS emits a quiescent thermal lu- 
minosity (BBR98) 
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where (M) is the time-averaged (over the NS core thermal 
timescale) mass-accretion rate onto the NS, and Q is the amount 
of heat depo sited in the crust per accreted nucleon (H aensel & 
Zdunik 1990; see pildsten & Rutledge 2000| for a review). For 
an accretion flux onto the NS (f accretion = eMc 2 /{4ttD 2 )), the 
"rock bottom" quiescent flux due to deep crustal heating is then 



Q 0.2 
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(2) 



ficiency e = 0.2 for accretion luminosity of (GM/R)(M). 

Since (M) = M out burst(^outburst/ T rec ), where T outbu i-st is the mean 
outburst duration and r rec is the mean recurrence timescale, this 
scenario relates the quiescent luminosity to the outburst prop- 
erties; comparisons to observations of several quiescent neu- 
tron stars (qNSs) match the predictions reasonably well (B rown 
et al. 1998; |Rutledge et al. 2000| ). It also helps us to understand 
why the black hol e systems in quiesc ence are so much fainter 
than NS systems ( Garcia et_al. 2001 ), as they cannot be ther- 
mally emitting. Equation (Q) is the minimum quiescent flux 
from t he neutron star. Sho u ld residual accretion occur in quies- 
cence ( fvlenou et al. 1999 ; Narayan et al. 1997), the accretion 
luminosity would be in addition to the thermal emission already 
present. 

KS 1731-260 is a transient, type-I X-ray bursting neutron 
star approxi mately 3.8 degrees from the Galactic center (V an 
Paradijs 1995). It was discovered in outburst in 1989, when its 
luminosity was 1.3 xlO 37 (d/8 kpc) 2 ergs" 1 ; subsequent anal- 
ysis found that it had been in outbur st at least since Oct 1988 
( Sunyaev 1989t jSunyaev et al. 1990 ), and has remained bright 
since then (see |Wijnands et al. 2001a for a complete description 
of observations). Sunyaev et al. (1990) suggested that it was 
a transient, apparently because it had not previously been de- 
tected, although no strong upper limits on previous X-ray emis- 



where (F mca - elion ) = t(M)c 2 /(4irD 2 ) is the accretion flux aver- 
aged over the NS core thermal timescale, and the accretion ef- 
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sion were given. (This was noted later by Barret et al. (1998) , 
who suggested the object may not be a transient at all. We 
adopt an outburst timescale of r ou tbuist ~ 13 yr, although the 
outburst may have been longer). Type I X-ray bursts were ob- 
served, at the rate of — 10 per day, and studied in detail with 
RXTE (|Muno et al. 2000|). Nearly coherent oscillations (580 



Hz) have been observed during the type-I bursts, and the fre- 
quency is interpreted as the spin frequency of the NS (S mith 
et al. 1997). The uncertainty in the localization of this source 
was initially 4.2' (90% confidence). An i mproved error circle 
(10" with ROSAT/HR1; |Barret et al. 1998[ ) found 13 candidate 
IR counterparts down to / = 15.4 (10 a). The integrated HI 
column density in the direction of KS 1731-260 is A^h,22=0.35 



(N H =N U .?.?. 10 22 cm" 2 ). iBarret et al. (1998)) found JV H ,22=1 



6 (A v =5-34; [Predehl & Schmitt 1995[ ), suggesting that N H is 
variable, and thus there must be significant contribution to A^h 
from the system itself. 

An RXTE/PCA scan of the Galactic Center found that 
KS 1731-2 60 had entered a low luminosity state (W ijnands 
et al. 2001a) in early 2001. [Wijnands et al. (2001a)| then ob- 
served the source with Chandra/ ACIS-S in late March (— 1 
months after the end of the outburst) and detected it at an unab- 
solved flux of- 2x lO-^ergcnrV 1 (0.5-10keV). 

1.1. Crust Thermal Emission 

The scenario described by equation (^) applies when the neu- 
tron star crust and atmosphere resemble that of a cooling neu- 
tron star - that is, the temperature decreases with increasing 
radius. For NS transients with short outbursts, such as Aql X-l 
(^outburst ~ 30 d), this is a good approximation, as the increase 
in cru st's temperature from the heati ng during the outburst is 
small (Ushomirsky & Rutledge 2001). This is not the case for 



KS 1731-260. The duration of its most recent outburst is of the 
order of (or longer than) the thermal diffusion timescale in the 
crust (BBR98). In a sense, KS 1731-260 in quiescence resem- 
bles a neutron star that accreted steadily at the outburst accre- 
tion rate, except that the core temperature will be at the lower 
value set by the time-averaged accretion rate (over the previous 
10 4 yr). 

To illustrate how the crust is so dramatically heated dur- 
ing a long outburst, consider the rise in temperature if no 
heat were conducted a way from the crust during the outburst. 
As we describe in § 3.3, the total heat capacity of the re- 



gion of the crust where most of the heat is deposited is C — 
5 x 10 35 erg K _1 . During the outburst, the total amount of 
heat deposited is QM outburst /m p w 10 44 erg and so, if no heat is 
conducted away from the heating region, the temperature there 
can rise to 2 x 10 8 K during the 13 year outburst. Even when 
thermal conduction is taken into account, the rise in crust tem- 
perature can still be > 5xl0 7 K, which is the typical temper- 
ature the crust and core would have if KS 1731-260 accreted 



steadily at (AT) = M ou tbur S t(T utburst/T re c) 
(for T lec Rj 1500 yr; see § 3.2). 



2.6 x 10- n M Q yr 



If the crust, which is composed of the ashes of H/He burn- 
ing, has a low thermal conductivity, then there is a substantial 
temperature gr adient between crust and core. As discussed by 



Brown (2000) for the case of steadily accreting neutron stars, 



when there is a substantial thermal gradient in the inner crust, 
the temperature of the crust reaction layers becomes decoupled 
from that of the core. 

W3nH, at ittp://heasarc.gsfc.nasa.go\ ; [Dickey & Lockman 199c| 
7 http://asc. harvard.edu/ciao2. 1/ 



For the NS in KS 1731-260, this means that, until the crust 
has thermally relaxed, we cannot directly infer its core tem- 
perature, as is possible for short-outburst transients such as 
Aql X-l. However, this also means that one can directly com- 
pute the current quiescent luminosity from KS 1731-260 using 
the fluence during the last outburst, as the thermal state of the 
crust is only weakly sensitive to the temperature of the core, 
and, hence, is nearly independent of the uncertainty in the ac- 
cretion history over the past — 10 4 yrs (the thermal time of the 
core), unlike in the case of short-outburst transients. Our sim- 
ulations, presented here, predict a quiescent luminosity which 
agrees with the observed value to within observational and the- 
oretical uncertainties (factor of a few) for the fiducial value of Q 
and outburst fluence. Moreover, this dependence of the thermal 
flux in the crust-dominated regime on M out burst, and not on (M), 
allows us to predict the evolution of L q , which can be directly 
confronted with observations. 

The domination of the quiescent luminosity by the cooling 
crust sets KS 1731-260 apart from the other neutron star tran- 
sients with short duration outbursts and makes it an ideal lab- 
oratory for separately measuring the thermal properties of the 
crust. 

1.2. Outline of the Paper 

In th is paper we report our analysis of the Chandra obser- 
vation ( [vVijnands et al. 200 la| ) of KS 1731-260 in quiescence 
and describe our simulations of the thermal state of its crust 
and core. We begin in § |2] by describing the Chandra obser- 
vation and spectral analysis of KS 1731-260. We include a 
description of possible IR counterparts (§ ^2j). After describ- 
ing the observations, we then discuss the implications for the 
crust and core thermal structure as outlined above. We first ap- 
ply (§ 3.2) the scenario in which the crust heating for any single 
outburst is negligible, i.e., we treat KS 1731-260 as if it were 
a short-duration outburst , such as Aql X-l. This analysis has 
been done prev iously by Wijnands et al. (2001a). We then ex- 
plain, in § 3.3, why this analysis is inapplicable to this source, 



and how the quiescent luminosity is determined by the ther- 
mal properties of the crust (its thermal diffusion timescale and 
heat capacity) for the next 1-30 yr. We present simulations of 
the quiescent lightcurve and demonstrate that monitoring obser- 
vations can const rain these properties, s imilar to proposals for 



glitching pulsars (Van Riper et al. 1991 ; Chong & Cheng 1994 
Hirano et al. 1997t Cheng et al. 1998[ ). We conclude in § p| with 



a summary and discussion of these results. 

2. Chandra OBSERVATIONS 

The data were obtained from the Chandra public archive. 
Th e discovery and observa tions have been previously analyzed 
by jWijnands et al. (2001a) , and details of the trigger and history 
of KS 1731-260 are included there. The observation was trig- 
gered when an RXTE/PCA scan of the Galactic Center found 
that KS 1731-26 had entered a low luminosity state (W ij- 
nands et al. 2001a), and was made 2001 Mar 27 0:18-06:23 
(TT) with the ACIS-S instrument ([Weisskopf 1988[) for a total 



exposure time of 19401 seconds. The X-ray source was imaged 
on the S3 (backside illuminated) chip, which was operated as 
a 1/4 array with 0.8 sec exposures. We analyzed the data us- 
ing the CI AO v2.1 7 with CALDB v2.6 and XSPEC vll (A r- 
naud 1996). 
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Two X-ray point sources separated by 31.3±0.1" are found 
with celldetect above a signal/noise ratio of 5, listed in Table [j]. 
We compared this image with an archived ROSAT/HRl image 
of this region (rh400718) 8 taken in 1997 when KS 1731-260 
was in outburst. The position of KS 1731-260 h ad been pre- 
viously determined this way (Barret et al. 1998) and clearly 
KS 1731-260 is at the position of Chandra source #1. 

X-ray source counts were extracted within an area 5 pixels 
in radius about the KS 1731-260 source position, for a total of 
183 counts. At 0.0075 counts/frame, the pileup fraction is neg- 
ligible. Background was taken from an annulus centered on the 
source position, with radii of 10 and 80 pixels. The expected 
num ber of background counts in the source region is 3. A KS 
test ( foess et al. 1995 ) finds that the times of arrival (TOA) of 
the 183 counts is consistent with a constant countrate. We place 
3cr upper-limits on variability of 3 x 1./V183 counts=22% rms. 
We also produced a power-density spectrum to search for a 
pulsed signal (number of frequency bins: 12750, with fre- 
quency resol ution of 5x 10~ 5 Hz, and a Nyquist frequency of 
0.625 Hz; see Press et al. 1995 ), using barycentered TOAs (tool 
axbary). No evidence for a pulsed signal is f ound: the largest 
Leahy-normalized power ( Leahy et al. 1983 ) was 20.37 (with 
a probability of chance occurrence from a Poisson-distributed 
countrate of 0.48). The absence of a coherent signal is not 
surprising in l ight of the detecti on of pulsations during type-I 
X-ray bursts ( |Muno et al. 2000| ) at ~580 Hz, well above our 
Nyquist frequency. 

2.1. Spectral Analysis ofKS 1 731—260 

We binned the KS 1731-260 data into 10 PI bins (0.5-10.0 
keV), and fit several single component spectral models (power- 
law, H atmosphere, or Raymond-Smith, a multicolor disk and 
blackbody). Galactic absorption is initially left as a free param- 
eter. The best-fit models were all statistically acceptable; the 
parameters are given in Table ||. 

While all models we investigated are statistically acceptable, 
some can be argued against on physical grounds. The pure 
power-law spectrum is unusually steep (a=5.2±0.6), typical 
more of a thermal spectrum than other processes. The emis- 
sion measure J n e ni,dV from the Raymond-Smith spectrum 
is higher by 2 orders of magnitude than typical from active 
stellar coron ae in the analogous RS CVn systems (D empsey 
et al. 199 3a; pempsey et al. 1993b|; see Bildsten & Rut- 
ledge 2000 for discussion on coronal emission from compan- 
ions in X-ray transients). The disk black-body spectral model 
implies an inner disk radius of ^0.7 km - considerably smaller 
than a neutron star. Ou r best-fit blackbody spe ctrum is consis- 
tent with that found by Wijnands et al. (2001a). 

As pointed out by BBR98, for accretion rates <2xl0~ 13 
M Q yr _1 , gravity stratifies metals in the NS atmosphere faster 
than they can be pr ovided by accretion (Bildsten, Salpeter, 
& Wasserman 1992), making a pure H atmosphere the ap- 
propriate description of the NS photosphere in quiescence. 
Due to the strongly energy-dependent opacity of free-free tran- 
sitions, these spectra are significantly different from black- 



bodies (Rajagopal & Romani 1996; Zavlin et al. 1996). Un 



like the results from black-body fits, the inferred NS radii 
(Roo=r / \ - 2GM / {rc 2 ) where r is the proper radius, and 
oo means the value as measured by a distant observer) of 
the four observed field transients using H atmosphere spectra 

8 http://heas arc.gsfc.nasa.gov 

'http://www.ipac.caltech.edu/2mass/overview/about2mass.html 



are consistent w ith what are expected t heoretically: Roo^ll 
km (Cen X-4, |Rutledge et al. 2001c|; Aql X-l, Rutledge 
et al. 2001a; 1608-522 [Ru tledge et al. 1999| ; and 4U~ 2129+47 
Rutledge et al. 200C), thus supporting the notion that much of 
the emission originates from a neutron star surface. 

The quiescent spectrum for qNSs is thus interpreted as ther- 
mal emission from a pure Hydrogen atmosphere NS, possibly 
with an underlying power-law (whose origin is not understood) 
which dominate s the spectrum at high (>3 keV) energies (see 
Rutledge et al. [200 lcf |2001a[ ). The spectral fit for a H at- 
mosphere spectrum alone (no power-law) gives &rjj?=120±30 
eV, an emission area radius of /? oo =10^5° (d/8 kpc) km, and 
an unabsorbed flux of 1.8x 10~ 13 ergcm" 2 s _1 (0.5-10 keV; the 
absorbed flux is 3.8xl0~ 14 ergcm" 2 s _1 ), corresponding to a 
luminosity of L x = 1.4xl0 33 (d/ 8kpc) 2 ergs-'(0.5-10 keV). 
The 90% confidence range in the unabsorbed flux is (1.1- 
4.6)xl0~ 13 ergcm" 2 s _1 . The best-fit spectrum is shown in 
Fig. [l]. The bolometric thermal luminosity (as observed 
at infinity) is L holoo = 4 7 rfl 00 2 cr(7;f ) 4 =2.7 xlO 33 (d/8 kpc) 2 
erg s" 1 (uncertain by a factor of ^3, due largely to spectral 
uncertainty). If 100% of the emergent luminosity were from 
accretion (at an efficiency of 0.2), this sets M = 2.3xl0~ 13 
M yr _1 , sufficiently low for the assumption of an H atmo- 
sphere. 

When we include an underlying power-law with all parame- 
ters free, the S/N of the spectrum does not constrain the model 
parameters to better than an order of magnitude. When we 
hold Roc and a fixed at typical values (/?oq=12.5 km for d= 8 
kpc; a = 0.85; cf. Rutledge et al. 2001cj ; |Rutledge et al. 2001a| ), 
the spectral model provides an acceptable fit to the data. This 
demonstrates that the spectrum of KS 1731-260 in quiescence 
is consistent with that observed from other qNSs. 

2.2. Possible Infrared Counterparts 

A comparison between Chandra astrometry and 2MASS im- 
ages finds a possible IR counterpart at the position of source #2 
(see below) but not of KS 173 1-260, with limits in J, H and K s 
of 15.8, 15.1 and 14.3 magnitudes respectively 9 . Comparing 
with IR sources previously examin ed as possible counterparts 
to KS 1731-260 ( [Barret et al. 1998[ ), the closest object is source 
"H", which is 1.5±1.1" away (we adopt 1" positional uncer- 
tainty for Chandra, and 0.5" uncertainty for the IR source). The 
second closest is "G", which is 3.3" away. In addition, there 
is a fainter IR object present in the H-band image ("new" ob- 
ject), but not in the J-band image taken by Barret, which is not 
marked, but which is consistent with the Chandra position for 
KS 1731-260. Further work in the IR - for example, searching 
for ellipsoidal variations in quiescence - is needed for a positive 
identification of the counterpart. 

Source #2 is consistent in position with 2MASS J173412.7- 
260548 (r = 0.46" in distance, about the accuracy of the 
2MASS position reconstruction error). For an IR field source 
density of p=6.9xl0" 3 arcsec" 2 (217700 objects in the 2M ASS 
catalog within a 1000" radius) the probability of a chance co- 
incidence is prob=l-exp(-47rr 2 p)=3.9%, which is the signif- 
icance of the association of the X-ray and IR objects. The 
source's 2MASS magnitudes and colors are J=12.264±0.033, 
7-i/=0.948 and7-/: s =1.33. 

When th is work was largely complete, Wijnands 
et al. (2001b) used relative astrometry with 2MASS, and found 
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that the Chandra localization is coincident with the "new" ob- 
ject and not "H" and "G". 

3. CORE DOMINATED EMISSION, CRUST DOMINATED EMISSION, 
AND THE FUTURE LIGHTCURVE OF KS 1731-260 

Having described the observations and spectral analysis of 
KS 1731-260, we now turn our attention to the thermal prop- 
erties of the neutron star's crust and core. We begin by esti- 
mating the amount of mass deposited during the previous out- 
burst; we then describe what the thermal structure would be if 
the chang e in crust temperature during this outburst were small 
(BBR98; Wijnands et al. 2001a ). This is the case for transients 
with short-duration outbursts. Because of the length of the out- 
burst for this source, however, the crust is substantially heated; 
§ 33 describes how the observed emission is determined by 
the heating of the crust during the past outburst, and not by the 
core equilibrium temperature. Using the method of UR01, we 
simulate the evolution of the quiescent lightcurve for different 
regimes of the crust and core microphysics. 



There is an uncertainty of about ±0.5 dex in this value, mostly 
due to the uncertainty in the bolometric value of F q . This esti- 
mate is simil ar to the one arrived at independently by Wijnands 
et al. (2001a), although our estimated mean recurrence time is 
longer (1500 ±0.5 dex yr vs. 200 yr), since our time-averaged 
outburst flux is a fac tor of six greater than that assumed by W i- 
jnands et al. (2001a). 

Strictly speaking, this estimate of T rec is the recurrence time 
to an outburst of comparable fiuence. A month-long outburst 
would have just 1% of the fiuence of the previous outburst and 
would not affect the estimate of eq. ([}]). This estimate also ne- 
glects any neutrino emission, which reduces the effective value 
of Q and hence r rec . More seriously, however, this analysis ig- 
nores the fact that the crust is strongly perturbed away from the 
core equilibrium temperature during such a long outburst. As 
a result, the quiescent emission currently observed is set by the 
thermal relaxation of the heated crust, and as we now demon- 
strate, is mostly decoupled from the thermal state of the core. 



3.1. The Outburst Fiuence 

KS 1731-260 has been observed several times, always in 
outburst, since 1988, and nearly continuously so since Jan 
1996 with the RXTE/ASM 10 . To estimate (F outburst ), we inte- 
grated the ASM counts using only data in which there was 
a > 3a 1-day detection and taking only 1 day detections in 
which the previous observation was made <4 days (to insure 
adequate sampling) prior. Periods in which the time to the 
previous observation was >4 days were not included. The 
total fiuence is 13660 ASM c/sxdays, over an integration of 
1599 days (including both detections and non-detections). If 
we assume an additional countrate at the 2a level on those 
observations when the countrate was below 3a detection, the 
total fiuence would only increase by 7%, which is our 2a 
upper-limit for the uncertainty in integrated counts. Assum- 
ing a £7b remss =5 keV bremsstrahlung spectrum with A^h,22=1-0 
(which gives 6.5xl0~ 10 erg cm" 2 s"'per 1 ASM c/s, corrected 
for absorption, 0.01-20 keV, from W3PIMMS 11 ; in the 2-10 
keV range, the conversion factor is 3. Ox 10~ 10 ergcm" 2 s"'per 1 
ASM c/s ), we find a mean outburst flux, corrected for absorp- 
tion, of 13660 ASM c/s x d/1599dx (6.5 xlO" 10 ergcm" 2 s" 1 
per 1 ASM count) = 5.6xl0~ 9 ergcm" 2 s" 1 . We conserva- 
tively estimate the spectral uncertainty to be at the 40% level 
(for a change in kT\, remss between 2 keV and 8 keV, the ASM 
counts/flux conversion changes by 40%). Extrapolating this 
mean to the entire outburst, we find an outburst fiuence of T = 
5.6x 10" 9 ergcm" 2 s _1 x (13yr x 3.2xl0 7 s/yr) = 2.3ergcm" 2 . 
This estimate is greater by a factor of ^6 than that of W ijnands 
et al. (2001a), who estimated T using the minimum observed 
outburst flux and an outburst duration of 11.5 yr, rather than 
13 yr. 

3.2. A First Estimate: Core-Dominated Emission 

The quiescent bolometric thermal flux is F q = 3.5 x 
10~ 13 ergcm" 2 s" 1 . If the change in crust temperature during the 
outburst were small, then equation (||) would apply: the time- 
averaged accretion flux is (F ac cretion) = 3~ ' /t kc , so that 



Q 



0.2 



\35F q 1.45 MeV 



1500yr. 



(3) 



3.3. Crust Dominated Emission 

During an outburst, the heating from the reactions raises 
the temperature in the crust around the reaction layers. Af- 
ter the outburst ends, the crust thermally relaxes and the ther- 
mal profile comes to resemble a cooling neutron star, i.e., the 
temperature decreases with radius. For neutron star transients 
with "typical" outburst fiuences and recurrence times (that is, 
■^outburst ~ 1 month, with an outburst luminosity ~ 0.1 Edding- 
ton), the variations in the crust temperature and hence L q are 
small (BBR98; UR01). For example, UR01 found that L q var- 
ied by ~ few% for T rec = 1 yr and by ~ 30% for T rec = 30 yr. 
The variation is small because the amount of energy deposited 
in such short outbursts is not significant compared to the heat 
content of the crust, and so the crust is not heated substantially 
compared to its temperature when in thermal equilibrium with 
the core. 

As our order of magnitude estimate of the crustal temper- 
ature for KS 1731-260 in § 1.1 demonstrates, the situation 
is completely different in the case of KS 1731-260, with its 
Toutburst ~ 13 yr outburst. 

To ascertain the magnitude of the deviation of L q from 
the value predicted by Eq. (jl|), we performed several simu- 
lations of thermal relaxation of the NS crust in a transient 
with observational properties (r ou tbursti r rec, ^outburst) inferred for 
KS 173 1-260 using the methods and microphysics as described 
in UROl. In summary, we solve the non-relativistic heat equa- 
tion in the crust, from p » 10 s g cm" 3 to 1.5 x 10 14 g cm" 3 . 
Heating due to nuclear reactions in the deep crust is simulated 
by depositing energy at densities corresponding to th e nuclear 
transitions computed by Haensel & Zdunik (1990), with the 



amount of energy set by the instantaneous accretion rate. At the 
outer boundary, we use the flux-tem peratu re relation for a fully 
accreted crust from Potekhin et al. (1997 ). The core is taken to 



be isothermal (a good assumption for the timescales of interest 
because of the large thermal conductivity), and its temperature 
evolves according to the mismatch between its neutrino emis- 
sivity and the flux from the crust. We start our simulations with 
the temperature profile corresponding to persistent accretion at 
the rate (M) corresponding to the time average of M utburst over 



10 http://xte. mit.edu 

1 'http://heasarc.gsfc.nasa.gov 
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the recurrence time. We then evolve the model through sev- 
eral outburst/quiescence cycles, as necessary for the model to 
"forget" the initial conditions and reach a limit cycle. 

If the crust resembles that of a neutron star steadily accreting 
at Moutburst, then the estimate of the recurrence time (eq. [|^]) is 
inapplicable (except as a lower limit). This is because the crust 
temperature is decoupled from that of the core. A full survey 
of parameter space should include simulations over a range of 
T rec . This is beyond the scope of this initial paper; rather, we 
presume that r rec is given by eq. (Jfy in order to survey the influ- 
ence of the crust and core microphysics. 

In these simulations, we examine the two main uncertain- 
ties in NS microphysics, (1) the impurity fraction, and, hence, 
the conductivity of the crust and (2) the possible presence of 
"enhanced" core neutrino cooling mechanisms, 12 such as direct 
Urea or pion condensation. First, the composition of accreting 
NS crusts is set by the nuclear processing of produ cts of burn- 
ing in the upper atmosph ere. Published calculations ( [Sato 1979 



Haensel&Zdunikl99CD presume that burning proceeds to pure 



iron, and, hence, the crust has no i mpurities. Howeve r, recent 
work by S chatz and collaborators dSchatz et al. 1999]; S chatz 
et al. 2001) shows that burning products are a mix of elements 
substantially heavier than iron, with the impurity parameter 13 
Qi mp comparable to the square of the average nuclear charge 
(Z) 2 . We therefore set boun ds on the crustal conduct ivity by 
using electron-ion scattering ( Yakovlev & Urpin 198C ), which 
has the same form as electron-impurity scattering with gimp = 



Z , and electron-phonon scattering (Baiko & Yakovlev 1995), 



which is appropriate for a pure crystal. We refer to these two 
cases as "low K" and "high K", respectively. Second, the uncer- 
tainty in core cooling is covered by simulations with "standard 
cooling" and "enhanced cooling" (modified Urea appropriately 
suppressed by nucleon superfluidity, and neutrino emission as 
when the core is a pion condensate, respectively; see UR01 for 
details). 

To illustrate the deviation of the crustal temperature from the 
equilibrium core temperature and the subsequent thermal relax- 
ation of the crust, we show in Fig. |2]the time evolution of the 
crustal temperature in the low K, standard cooling case for ob- 
servational parameters inferred for KS 1731-260 (r rec = 1500 
yr, T outburst = 13 yr, and M ou tbur S t = 3 x 1O" 9 M yr" 1 ). Verti- 
cal slices of the surface in the density-temperature plane show 
the instantaneous temperature distribution in the crust, while 
slices in the time-temperature plane show the time evolution of 
the temperature since the beginning of the outburst. Clearly, 
the crust is heated to well above its pre-outburst temperature 
profile, which tracks the core temperature (note that, because 
the crust is heated to well above its temperature in equilib- 
rium with the core, this qualitative result is not very sensi- 
tive to the assumed r rec ). A steep temperature gradient carries 
most of the heat from the deep crustal heating region around 
10 12 -10 13 g cm" 3 into the core, while a shallower temperature 
gradient carries a fraction of the heat towards the surface. The 
temperature near the top of the c rust (^10 10 g cm" 3 ; a nd, hence, 



the X-ray luminosity, L q oc T^ 4 ; Potekhin et al. 1997| ) reaches a 



maximum at the end of the outburst (time=13 years), and then 
decays back to the pre-outburst value on the thermal timescale 
of the crust (~ 30 yrs for this model). Thermal evolution in 
other cases is similar (but, of course, is different in the abso- 



lute magnitude of the temperature change), and the temperature 
changes are always much bigger than those that occur during 
short-duration outbursts (cf. Fig. 1 of UR01). 

In Fig. ||, we show the long-term evolution of the quiescent 
surface luminosity from KS 1731-260 after accretion onto the 
compact object has ceased for the four cases discussed above, 
namely, low and high crustal conductivity K, and standard and 
enhanced core neutrino emission (see the legend on the plot). 
First consider Fig. ||a, where we set Q = 1.45MeV. In all four 
cases, the transition from crust-dominated (at early times) to 
core-dominated cooling (at late times) is evident as a drop in the 
luminosity, ranging from 50% (high K, standard cooling case) 
to a factor of 100 (low K, enhanced cooling case). Regardless 
of the core neutrino emissivity, this transition occurs at ~ 30 yrs 
for low conductivity crusts, and after ~ 1 yr for high conductiv- 
ity crusts. These timescales are easy to understand, as they are 
just the thermal time to the appropriate depth in the crust (cf. 
Figure 3 of UR01). Since the low K crust has a longer thermal 
time, it stays hot about 30 times longer than the high K crust. 

After the crust thermally relaxes, L q is set by the emission 
from a hot core. If the only core neutrino emission mechanism 
is modified Urea, suppressed by nucleon superfluidity, then the 
amount of heat lost by neutrinos from the core is negligible. 
Therefore, all the heat deposited into the star comes out as 
thermal emission, and L q asymptotically approaches the value 
given Eq. ([]]), with Q w 1 .45MeV (see solid and dashed lines 
in Fig. g). On the other hand, when enhanced neutrino cooling 
i s allowed, ~ 90% of the deposited heat escapes as neutrinos 
( |Colpi et al. 2000j ; UR01), and only the remaining ~ 10% is ra- 
diated thermally from the surface. In this case, the "effective" 
value of Q in Eq. ([!]), in the sense of the energy retained in the 
star and re-radiated thermally, is Q w O.lMeV. 

As shown above, the crust is heated to well above the core 
temperature during the long outburst. How does the tempera- 
ture difference between the crust and the core depend on the 
crustal conductivity? As an estimate, we can assume that, dur- 
ing the outburst, the crust comes to an equilibrium where the 
heat deposited around p rs 10 12 - 10 13 g cm" 3 at a rate E is con- 
ducted away from this region. This is a very good assumption 
for the high conductivity case, and an acceptable one for the low 
conductivity case, where the crustal thermal time is comparable 
to the outburst duration. During the outburst, most of the heat 
is conducted into the core with a flux F sa K(T cm $t - T core )/ Ar, 
where Ar w 500 m is the distance between the crustal heating 
region and the core. Thermal balance then requires E/4irR 2 sa 

^(^crust ^core)/Ar, Or 7cmst ^core ~ 

1.2 x 10 8 K /K 19 , where 
Kig is the conductivity in units of 10 19 erg cm" 1 s K" 1 . In the 
low conductivity case, K\g ~ 1, so r crust 3> r c0 re- We see from 
Fig. ||a that L q at early times in the low K case is the same for 
standard and enhanced cooling (solid and dotted lines). This is 
now easy to understand, since, as we showed above, the crustal 
temperature is decoupled from the core temperature if the con- 
ductivity is low. This situation is ex actly t he same as for persis- 
tent accreters discussed by Brown (2000). On the other hand, 
for high conductivity crusts, ^19 > 10, so r clust can not signifi- 
cantly deviate from T COK during the outburst. Therefore, crust- 
dominated L q in the high K case (dashed and dash-dotted lines 
in Fig. |^a) more closely reflects the core temperature, which 
is quite different between the standard and enhanced cooling 



12 Neutrino emission due to crustal bremsstrahlung is insignificant at the relevant temperatures, but is included in all models. 
13 The impurity parameter is the deviation of the nuclear charges Z of the elements of the mix from the average charge (Z), 2imp : 
takes into account the relative abundance of the species. 



((Z-(Z)) 2 ), where the average 
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From examination of lightcurves in Figs. ||a, it is likely that 
KS 1731-260 was observed in quiescence during the crust- 
dominated phase, and is presently evolving through this phase. 
In three of four cases considered in Fig. ga, the predicted qui- 
escent luminosity L q right after the outburst is a factor of 3 
higher than the observed L\, Q \ (itself uncertain by a factor of 
3, due to spectral uncertainty). Since L q depends directly on 
Q J Mouibm-stdt over the outburst duration, this difference could 
also be due to a systematic overestimate of M ou tburst during the 
« 13 yr outburst, only the last 5 years of which have been cov- 
ered ~ daily with RXTE/ASM. Alternatively, this difference 
could be due to a value of Q different from the fiducial, since 
the exact value of Q is uncertain (see Schatz et al. 1999, 2001 
for a discussion on different crustal compositions than assumed 
by Haensel & Zdunik 1990| ). To provide lightcurves which ex- 
trapolate from the currently observed luminosity, in Fig. ||b we 
adjusted Q (or, equivalently, / Moutburst ^0 by factors of 0.3 to 2, 
as noted in the figure caption, to obtain L q ~ 2.7 x 10 33 erg s -1 
at the start of quiescence. 

We stress that the overall normalization of these lightcurves 
is uncertain due to uncertainty in Q and outburst fluence, but 
the shape is determined by the microphysics of the crust and the 
temperature of the core. The time evolution of the luminosity 
will permit distinction between the four cases, and, hence, di- 
rectly infer the integrated conductivity of the crust and possibly 
indicate whether enhanced neutrino emission mechanisms are 
operating in the core. In particular, observing the drop in the lu- 
minosity will permit the measurement of the thermal timescale 
of the crust. In addition, the relative magnitude of the drop will 
tell us about the presence or absence of enhanced neutrino emis- 
sion from the core. Such measurements (accurate to 510%) are 
well within the capabilities of present X-ray instrumentation, as 
the present detection demonstrates. Observations over the next 
year will be able to exclude the case where the crust is a pure 
crystal and hence has very high conductivity. 

Since the initial understanding of Type I bursts as thermonu- 
clear flashes, there has been an open question as to how much 
flux is ris ing up from deep parts of the sta r into the burning re- 
gion (see Paczynski 1983; feildsten 1995 ). Large values of this 
flux (~ 1 MeV per accreted nucleon) can sta bilize burning a t ac- 



cretion rates lower than otherwise expected. Brown (2000) first 
showed that most of the crustal heating in a persistent source 
went into the NS core and that, at most, 100 keV per accreted 
nucleon came out into the upper parts of the star. The quiescent 
observations of KS 1731-260 shortly after outburst strongly 
support this picture and tell us that deep nuclear heating will 
not impact the Type I X-ray burst behavior. 

Assuming that the source has been transiently accreting for 
the past 10 4 yr or so, with outburst fluences similar to that of the 
most recent observed outburst, we can constrain the core tem- 
perature to be less than the maximum crust temperature. To de- 
termine the most conservative (i.e., largest) upper bound on this 
peak crust temperature, we use the following assumptions. We 
take the atmosphere to be composed of pure iron (which has a 
higher opacity than hydrogen/helium) and we use low thermal 
conductivity (electron-ion scattering) in the crust. These as- 
sumptions lead to the largest temperature rise between the pho- 
tosphere and the crust. Taking kT e fs = 130eV, we find a peak 
crust temperature, and hence a maximum core temperature, of 
3.5xl0 8 K. 



4. DISCUSSION AND CONCLUSIONS 

The large amount of heat deposited during the recent (and 
the only observed) outburst of KS 1731-260, and the lengthy 
duration of the outburst compared to the thermal diffusion time 
of the NS crust, imply that the quiescent flux recently observed 
is dominated by emission from the cooling crust and not the 
core. The importance of this result is that the quiescent lumi- 
nosity is calculated using the properties of only the most recent 
outburst, rather than by estimating the accretion history over 
the past 10 4 yr, as is the case i n the core-dominated transients 
such as Aql X-l and Cen X-4 (|Rutledge et al. 200 lj; Rutledge 
et al. 2001a). Our prediction of the quiescent luminosity from 
the cooling crust following the outburst agrees, to within obser- 
vational and theoretical uncertainties (a factor of a few), with 
the estimate of the observed bolometric luminosity. Unlike the 
case of the short r ou tburst transients where we can directly infer 
the average core temperature, in the case of KS 1731-260 we 
cannot measure the core temperature, and hence the recurrence 
time, until after the crust has thermally relaxed (1-30 yr after 
the end of the outburst). Assuming that the source has been 
transiently accreting for the past 10 4 yr or so, with outburst flu- 
ences comparable to that of the most recent outburst, we can 
constrain the core temperature to be les s than the peak crust 
temperature, which, using the relation of Potekhin et al. (1997) 
for an iron crust (low conductivity), gives r core < 2x 10 B K. 

We predict future lightcurves for KS 1731-260 on the ba- 
sis of 4 different scenario for crustal and core microphysics 
(high/low conductivity, and enhanced/normal neutrino emissiv- 
ity). These can be distinguished from each other observation- 
ally - one in < 1 yr (excluding or confirming the high K crust 
with enhanced neutrino emissivity core), others over longer 
terms (10-100 yrs). This offers the unique opportunity to di- 
rectly constrain the properties of the NS crust - thermal con- 
ductivity, heat capacity, and depth to the NS core - through 
monitoring observations that would measure the predicted de- 
caying lightcurve. 

KS 1731-260 is now the fifth— after Cen X-4, Aql X-l, 
4U 1608-522, and 4U 2129+47— field transient neutron star 
which displays in quiescence a spectrum consistent with a 
H atmosphere NS of radius Roo^ 12 km, and (possibly) a 
hard (a ~0.85) power-law component which comprises 10- 
40% of the quiescent (0.5-10 keV) flux. KS 1731-260 shows 
no variability on short (<hours) timescales. These spectral 
and variability properties are identical to those of the well- 
studied Aql X-l and Cen X-4 dRutledge et al. 200 lcj; Rutledge 
et al. 2001a), and the spectral properties are similar to those 
from lower S/N data of 4U 1608-522 and 4U 2129+47 (R ut- 
ledge et al. 1999; [Rutledge et al. 2000f The derived radius is 
also comparable to that derived from this object using the in- 
dependent jnedio^_^f_a_£o^hn^adms_£xrja^ 
et al. 1997). 

It is possible that a persistent source can turn off after long 
periods of accretion (here, a persistent source is one accret- 
ing for a period ^10 4 yr, the core heating timescale). While 
Eq. 1 would seem to be applicable to predict the post turn- 
off quiescent thermal luminosity from the crust, the equation 
assumes that all heat deposited in the crust is re-radiated dur- 
ing quiescence. That will not be the case for sources ap- 
proaching Eddington luminosity, such as KS 1731-260, where 
the effect of neutrino cooling in the crust and core cannot 
be neglected. In such cases, L q (2crust {M omhmstl nip) = 9 x 
lO 33 (e crust /15OeV)(M outburst /lO- 9 M yr- 1 ). Here g cmst is the 
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energy per accreted baryon emitted through the atmosphere 
from the crust immediately following the outburst. A minimum 
value, 2cmst = 0.1(2 w 150 eV/nucleon, applies when the source 
accretes persistently at the Eddington limit, in w hich only 10% 
of the deposited heat flows towards the surface (Brown 2000). 
At lower persistent luminosities, however, the crust and core 
temperatures are lower and neutrino cooling in each becomes 
less important, so Q cmst increases. In the case of KS 1731-260 
(^outburst ~ O.lLEdd), we obtain Q crust -250 eV, if KS 1731-260 
were a persistent accreter. However, the quiescent luminosity of 
KS 1731-260 gives instead Qcmst = 20 eV, a factor of 7 below 
the lower limit, even for a persistent source accreting at the Ed- 
dington rate. This demonstrates that KS 1731-260 cannot be a 
persistent accreter, because its crust (and by extension the core) 
is too cool to have been accreting for ^10 4 yr. Note, however, 
that we have neglected the effect of of neu trino emission due to 



Cooper pairing of superfluid particles (see Yakovlev et al. 2001 



for a review) whic h may reduce the value of Qcmst- 

As discussed by Rutledge et al. (2000), the study of qNSs in 
globular clusters has yielded a new way to accurately measure 
NS radii. The distances to some gl obular clusters have been de- 
termined to <2% post-Hipparcos ( |Carretta et al. 2000|), and can 
be measured even more directly with SIM (Unwin 1998), ef- 
fectively removing distance as an uncertainty. At present, there 
is only one globular cluster qNS in which the thermal com- 
ponent has been s pectroscopically analyzed (in uj Cen; R ut- 
ledge et al. 2001b) and th ere are two more proposed (in 47 
Tuc; |Grindlay et al. 2001 ). Measuring NS radii to an accu- 
racy of ±0.5 km — even without knowing the mass — can ex- 
clude ~ 50% of the proposed equa tions of state for the NS core 
matter (Lattimer & Prakash 2001). Even though we presently 



cannot make such accurate radius measurements in the field 
transients due to uncertain distances (which can be overcome 
with SIM), they remain important targets for detailed studies of 



the observational phenomena of qNSs, such as we report here. 

The relative luminosity of the (possible) power-law compo- 
nent poses a puzzle. It has been suggested that this compo- 
nent is due to accretion onto the NS magnetosphere (C ampana 
et al. 1998). However, such accretion would be unrelated to 
the quiescent luminosity of NSs as predicted from deep crustal 
heating (which is dominated by outburst accretion), and it can 
only be ascribed to coincidence that in the cases where adequate 
sensitivity is obtained, the power-law component produces a 
comparable fraction of the thermal flux in different sources. 
This suggests that the power-law component and H atmosphere 
spectral component are perhaps more closely related than pre- 
viously thought. 
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FIG. 1 . — The vFv model spectrum of KS 1731-260, and the observed Chandra/ ACIS-S BI data. The solid line is the best-fit unabsorbed spectrum (H atmosphere 
plus power-law, with i?oo=12.5 km and o=0.85 held fixed) . The dashed line is the H atmosphere component, the dashed-dotted line is the power-law component. 
The crosses are the observed Chandra data, with error bars in countrate. The two spectral components are equal near 3.5 keV, above which the power-law component 
dominates, and below which the H atmosphere component dominates. 



FIG. 2. — Crustal temperature distribution as a function of time, for r ol nb urst = 13yr, 3 X 10~ 9 Mq yr -1 , recurrence time of 1500 yr, and deep crustal heating of 
Q = 1 .45 MeV, using "low K" conductivity, and modified Urea neutrino emission from the core. The x axis is the density, y axis is time since the beginning of the 
outburst, and z axis is temperature. 



FIG. 3. — Left panel (a): time evolution of the quiescent luminosity assuming an outburst accretion rate of 3 X 10~ 9 Mq yr -1 , recurrence time of 1500 yr, and 
deep crustal heating of Q = 1.45 MeV/accreted baryon. Solid and dotted lines are for low crustal conductivity (set by electron-impurity scattering), dashed and dash- 
dotted lines are for high crustal conductivity (electron-phonon scattering). Solid and dashed line correspond to simulations with standard neutrino cooling processes 
(modified Urea suppressed by nucleon superfluidity), dotted and dash-dotted lines - enhanced neutrino emission due to a core n condensate (with superfluid 
correction factors). Right panel (b): same as (a), except the nuclear energy release Qnuc in the crust is adjusted to give L q ?s 2.7 X 10 33 erg s~' just following the 
outburst. g nllc values are: solid line — 0.42 MeV, dotted line — 0.49 MeV, dashed line — 0.6 MeV, dash-dotted line — 3.1 MeV. 
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Table 1 

Detected X-ray Sources in the field of KS 1731-260 





R.A. 


Dec. 


Counts 


No. 


(J2000) 


(J2000) 




1 (KS 1731-260) 


17h34ml3.45s 


-26d05ml8.7s 


180±13 


2 


17h34ml2.68s 


-26d05m48.3s 


50±7 



NOTE. — Positional uncertainty (from systematics of spacecraft pointing; frittp://cxc. harvard.edu/mta/ASPECT/celmon; ) is 
Total background subtracted counts in 19401 sec integration. 



Table 2 

Chandra SPECTRAL MODEL PARAMETERS OF 

KS 1731-260 (0.5-10 keV) 



H Almospher 



"H.22 
M-°0(eV) 

«oo (km) 
Tolal Model Flux 
xj/dof(prob) 



120±30 



H Aim. + Power Law 



W H,22 
H3f(eV) 
ffoo (km) 



Tolal Model Flux 
xj/dof(prob) 



1.3 ±0.3 

90! ?a 

7 ,+30 



2.0 



H Aim. + Power Law Iff o 



™H.22 
ffoo (km) 



1.06±0.08 
111 ±3 
(12.5) 
(0.85) 
0.2 
2.1 



Pholon Power Law 



W H,22 

Tolal Model Flux 
xj/dof(prob) 



1.7±0.4 
5.2±0.6 
20.0 
0.94/7 (0.48) 



W H,22 
Z(Z sol ) 
kT (keV) 

J" n e nffdV 
Tolal Model Flux 
xj/dof(prob) 



0.7 



,+0.2 
-0.1 
(1) 
l.4±0.1 

(6.2±1.5)x 10 55 

0.93 

2.32/7 (0.02) 



Multicolor Disk 



%,22 
Tm (eV) 

ff in */cos(0) (km) 
Tolal Model Flux 
xj/dof(prob) 



1.1±0.2 
370±50 



%.22 
"ef? <eV) 
ffoo (km) 
Tolal Model Flux 
xj/dof(prob) 



0.90±0.2 
300±40 



Note. — X-ray fluxes are un-absorbed, in units of 10 erg em 
Values in parenthesis are held fixed. Assumed source distance d-8 kpc. 



s (0.5-10 keV). Uncertainties are lcr. 



